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ABSTRACT 


Nomographs  and  tables  based  on  the  Poisson  Distribution  are  provided.  The  graphs 
give  the  mean  time  between  faiures  (MTBF)  from  test  data  at  one  sided  confidence 
boundaries.  Operating  characteristic  (O-C)  curves  for  fixed  time  life  test  and  fixed 
sample  attribute  test  can  be  calculated. 
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LIST  OF  SYMBOLS 


Confidence  (I-Pa) 

Parameter  of  Poisson  Distribution 
Mean  Time  Between  Failure 
Number  of  Items  Tested 
Probability  of  Acceptance 
Fraction  Defective 
Number  of  Failures 
Test  Time 
Chi  Square 

(Sigma)  Summation  from  r^  to  r2 

Factorial  sign  01  =  1 
Reliability 


INTRODUCTION 


Components  in  missiles,  launchers  and  other  tactical  weapons  are  subject  to 
failure  during  their  useful  life.  It  is  frequently  required  to  calculate  confidence 
boundaries  from  life  test  data  or  from  attribute  tests. 

Two  methods  are  given  below  to  allow  rapid  and  easy  computation  of  these 
estimates  -  rearranged  Poisson  tables  and  nomographs.  The  tables  are  arranged  in  such 
a  way  as  to  give  the  Poisson  parameter  (a)  for  confidence  levels  ranging  from  2%  to 
99.9%  for  up  to  50  failures.  The  nomographs  encompass  two  distinct,  limited,  ranges  (2 
-30%  and  70  -  98%)  of  confidence  for  up  to  45  failures.  They  solve  the  equations  for 

_  T  or  a  =  NP 

O’  5  MTBF 

By  simply  using  a  straight  edge  connecting  line  A  (a)  with  B  (T  or  N)  and  reading  the 
answer  on  line  C  (MTBF,  P,  R). 

THEORY 

The  tables  are  based  on  the  Poisson  Equation 
P  «  51  are“a  1-C 

a  ito  7 — 

where  a  T  or  a  =  NP 

MTBF 

This  methodology  can  be  used  to  obtain  approximate  confidence  boundaries  for 
reliability  from  attribute  test  data  and  also  O-C  curves  by  using  the  Poisson 
approximation  for  the  binomial  equation,  letting 

a  _  T  NP  =  N(l-R)  *  3. 

"  MTBF 


It  should  be  pointed  out  that  Chi-square  tables  are  normally  used  to  calculate  the 

confidence  boundaries  by  the  equation 

MTBF  ^  2T 

X2  (2r  +  2,  1-C) 

Which  is  equivalent  to  the  Poisson  method  used  here. 

Since  the  Poisson  distribution  is  only  an  approximation  for  the  Binomial-not  very 
accurate  for  small  sample  sizes-a  nomograph  published  by  Western  Electric  Company 
could  be  used. 

f. 

1  See  appendix  for  a  description  of  the  computer  program  used  to  calculate  the 

t 


tables. 


USER'S  MANUAL 


The  following  examples  will  show  how  to  use  the  nomograph  and  the  tables: 

EXAMPLE  1: 

Problem:  A  sample  of  100  was  randomly  selected  from  a  population  whose 
reliability  we  wish  to  estimate.  There  were  three  failures.  It  is  desired  to  find  the 
lower  80%  confidence  limit. 

Solution:  Using  the  nomograph  (Figure  lb)  go  to  80%  confidence  on  the  baseline, 
follow  the  vertical  to  three  failures  -  go  parallel  to  the  baseline  over  to  the  reference 
line  A,  mark  this  point,  connect  it  with  100  on  line  (B)  and  extend  to  line  C  marked 
"fraction  defective"  -  read  0.055.  The  dark  line  is  the  worked  solution.  The  confidence 
limit  can  also  be  determined  using  the  Poisson  tables  and  the  formula. 

P  =  a/N 

where 

a  =  Poisson  parameter 
N  =  number  of  sample  units 
P  =  fraction  defective 

Using  Poisson  table  1,  sheet  2  read  5.5153  for  80%  confidence  and  3  failures. 

P  _a_  5.5153  0.055153 

N  100 

This  means  the  reliability  is  94.5%  with  80%  confidence.  The  "Binominal  Reliability 
Tables"*  give  .94554  for  the  reliability.  This  compares  favorably  with  our  0.05515. 

EXAMPLE  2: 

Problem:  An  equipment  should  have  a  Mean  Time  Between  Failure  (MTBF)  of 
10,000  hours,  but  is  acceptable  if  90%  have  a  MTBF  of  2000  hours  or  more  ((>  =  0.1). 
If  the  actual  MTBF  is  10,000  hours  the  producer  does  not  want  it  rejected  more  than 

♦Binominal  reliability  table,  NAVWEPS  Rep.  8090  (NOTS  TP  3140) 


10%  («t*  0.1)  of  the  times.  To  do  this  find  the  correct  test  time,  draw  an  O-C  curve, 
and  also  determine  what  the  O-C  curve  would  look  if  the  test  time  and  the  allowable 
failures  were  tripled. 

Solution;  Using  table  I,  sheets  1  and  2  find  the  test  time  for  10%  confidence  for 
10,000  hours  MTBF  and  90%  for  2000  Hours  MTBF.  (a  =  T/MTBF) 


Table  I,  Sheet  I  T  Table  I,  Sheet  3  T 


FAILURES 

10%  Conf 
a 

a  x  10,000 

90%  Conf 
a 

a  x  2000 

0 

.1056 

1056 

2.307 

4614 

1 

.5309 

5309 

3.8901 

7780 

2 

1.1020 

11020 

5.3229 

10646 

3 

1.7439 

17439 

6.6812 

13362 

4 

2.4319 

24319 

7.9937 

15588 

5 

3.1512 

31512 

9.2751 

18550 

6 

3.8946 

38946 

10.5325 

21064 

This  table  shows  that  the  test  time  most  closely  satisfying  both  conditions  is  11,000 
hours  with  two  failures. 

The  O-C  curves  for  2  and  6  failures  are  shown  on  the  next  page.  This  example  is 
illustrated  in  the  "Pocket  Handbook  on  Reliability,"*  Section  5,  pp.  52-63.  It  should  be 
noted  that  the  calculated  points  close  to  the  2000  and  10,000  hour  MTBF  lines  coincide 
with  the  values  given  in  Figure  5-4,  p.  58. 


•"Pocket  HandbooK  on  Reliability"  by  J.K.  Byers,  and  H.A.  Weibe,  USAAVSCOM  Tech 
Rep.  77-16. 


Use  the  nomograph  for  70,  80,  90  and  95%  confidence 


Pa 

Confidence  (2  failures) 

MTBF  at  11,000 

Tables 

.3 

70 

3050 

3042 

.2 

80 

2550 

2570 

.1 

90 

2050 

2066 

.05 

95 

1750 

1797 

Use  the  Tables: 

Pa 

Confidence  (2  failures) 

MTBF  at  11,000 

.9 

10 

9981 

.8 

20 

7161 

.7 

30 

5748 

Probability  of  Acceptance 


1000  2000  3000  4000  5000  6000  7000  8000  9000  10,000 

Mean  Time  Between  Failures 


EXAMPLE  3: 


The  "Pocket  Handbook  on  Reliability"  on  pages  40-42  illustrates  a  method  for 
calculating  the  MTBF  for  certain  confidence  limits.  It  gives  the  lower  95%  confidence 
limit  for  1836  hours  and  3226  hours  of  testing  with  4  and  7  failures,  respectively,  as 
200.6  and  245  hours.  Using  the  confidence,  the  test  time,  and  the  number  of  failures 
and  Nomograph  Figure  lc  gives  these  answers  by  going  from  95%  confidence,  4  failures 
parallel  to  the  baseline  to  line  A.  Connect  the  point  on  line  A  with  1836  on  line  B  for 
total  test  time,  extend  and  read  200  MTBF  on  line  C.  Use  the  same  method  for  3226 
hours  and  7  failures. 

EXAMPLE  4: 

Problem:  A  piece  of  equipment  was  tested  for  10  hours  and  had  3  failures.  What 
is  the  expected  MTBF  at  80%  confidence. 

Solution:  Go  to  80%  confidence  on  the  baseline  of  Nomograph  Figure  la  follow 
the  vertical  to  3  failures,  go  parallel  to  the  baseline  over  to  line  A,  mark  this  point, 
connect  with  10  on  line  B,  extend  to  line  C,  and  1.8  on  the  right  side  of  li,.e  C. 


APPENDIX 


A  computer  program  is  used  to  calculate  the  probability  based  upon  the  following 
equations: 


r=50 


r  -a 
a  e 

r! 


a*_I_ 

MTBF 
a  =  NP 


When  Pf  =  probability  of  success 
r  =  number  of  failures 


r  =  50 


sigma  summation  from  0  to  r 


a  =  Poisson  parameter 
T  =  total  test  time 
MTBF  =  Mean  Time  Between  Failure 
N  =  number  of  items  tested 
P  =  fraction  defective 

The  calculational  scheme  is  developed  making  the  following  expansion: 


0  -a  _1  -a  2  -a 


P_  a  e  _a'e  ”  ^a_e 

r  nr  ^ 


0! 


1! 


2! 


-  + 


r  -a 
a  e 


removing  the  common  term  e  >  the  equation  then  becomes 


•  •  • 


the  equation  becomes 

Pr  =  e”a(DQ  +  Dj  +  D2  +  •  •  •  +  Dr_j  +  Dr) 


With  this  equation  a  table  of  probabilities  for  a  from  0.01  to  100  in  increments  of  0.1 
0.1  and  for  failures  from  0  to  50  in  increments  of  1  is  created  as  shown  in  Table  A. 

It  is  then  possible  to  enter  Table  A  with  a  failure  and  a  probability  and  interpolate 
for  a.  The  results  are  shown  in  Table  B.  It  is  from  the  use  of  Table  B  that  the 
nomographs  are  created.  For  direct  sampling  a  is  the  product  of  the  number  of  items  in 
the  sample  and  the  fraction  of  failures  (i.e.,  a  =  NP).  For  time  testing  a  is  the  test  time 
divided  by  the  MTBF  (i.e.,  a  =  T/MTBF). 

The  nomographs  do  not  show  probability  but  "confidence."  Confidence  is  1.0 
minus  the  probability  of  failure.  This  means  if  the  probability  of  failure  is  only  .1  that 
the  MTBF  will  be  lower  than  the  specified  amount  then  the  confidence  will  be  .9  that 
the  MTBF  will  be  higher  than  the  specified  amount. 

The  points  on  the  left  side  of  the  nomograph  give  the  numerical  values  of  a.  This 
value  then  is  used  to  calculate  the  total  test  time  (T)  or  MTBF  in  the  formut* 

T  (or  the  fraction  defective  in  the  formula  a  =  NP). 

MTBF 

Since  the  Poisson  distribution  is  only  an  approximation  for  the  binomial  -  not  very 
accurate  at  small  sample  sizes  -  a  nomograph  of  the  binomial  published  by  Western 
Electric  Company  could  be  used.  Also,  Sandia  Corporation,  in  April  1960,  published  "A 
Confidence  Limit  Computer"  based  on  the  binomial  distribution. 

Table  C  shows  the  input  data.  The  computer  program  is  attached  as  Table  1. 


PROBABILITY 


1  2  ...  48  49  50 


0.9048  0.9953 

0.8187  0.9825 

0.7408  0.9631 

0.9998 

O.80/8 

0.9964 

1.000 

1.000 

1.000 

1.000  1.000 

1.000  1.000 

1.000  1.000 

Probabilities 

TABLE  A 
for  Lambda 

Versus  Failures 

0.98 

0.95 

0.05 

0.02 

0.0127 

0.0476 

2.9959 

3.9126 

0.2128 

0.3531 

4.7448 

5.8349 

36.5707 

38.9644 

62.1713 

65.5711 

37.4261 

39.8484 

63.2872 

66.7158 

TABLE  B 

Lambdas  for  Failures  Versus  Probabilities 


SAMPLE  INPUT 


0.95  0.90  0.S0 

0.60  0.50  0.40  0.50 

0.10  0.05  0.02 

0.005  0.001  0.0005  0.0001  0.00005  0.00001 


prob 

.98  I'OO 

.95000 

.90000 

.80000 

SCONF 

2.00  000 

5.00000 

10.00000 

28.00000 

ALPHA 

.  02  000 

.05000 

.10000 

.20000 

.0127 

.0476 

.1056 

.2240 

.2128 

.3531 

.5309 

.8242 

.5643 

.8167 

1.1020 

1 .5347 

1.01 50 

1.3650 

1.7439 

2.2967 

1.52  79 

1.9690 

2.4319 

3.0894 

2.0884 

2.6125 

3.1512 

■•.9036 

2.6  8  32 

3.2848 

3.8946 

3.30  68 

3.9802 

4.6555 

5.5758 

3.95  18 

4.6951 

5.4320 

f  .4282 

4.61 77 

5.4248 

6.2209 

7.2891 

5.30  00 

6.1683 

7.0204 

°.1567 

5.9957 

6.9237 

7.8289 

a 

.0307 

6.7041 

7.6893 

8.6455 

o 

.9100 

7.4231 

8.5633 

9.4692 

in 

.7939 

8. 15  22 

9.2457 

1  0.2996 

1 1 

.6819 

8.8911 

10.0354 

11.1349 

12 

.5737 

9.6  3  69 

10.8316 

11.9758 

r 

.4689 

1 0. 39  12 

11.6338 

12.8214 

14 

.3673 

1 1.15 12 

12.4414 

13.6712 

1  K 

.2685 

1 1.91 83 

13.2541 

14.5250 

14 

.1723 

12.69  11 

14.0716 

15.3825 

17 

.0786 

13.4687 

14.8936 

16.2432 

1? 

.9871 

14.2516 

15.7192 

17.1075 

10 

.8977 

15.0373 

16.5486 

17.9743 

1  c 

.8102 

15.8314 

17.3818 

18.8440 

2n 

.7245 

16.62  76 

18,2183 

19.7165 

21 

.6405 

17.42  77 
18.2314 
19.0386 
19.84  91 
20.66  29 


21.4797 
22.29  96 
23.12  15 
23.9462 
24.7737 


19.0577 

19.9006 

20.7455 

21.5939 

22.4441 


23.2974 

24.1523 

25.0100 

25.8693 

26.7308 


20.5914 

21.4684 

22.3477 

23.229? 

24.1127 


24.9981 

25.8851 

26.7738 

27.6642 

28.5562 


22.5582 
?  3  •  4  7  7  5 
24.3983 
?5 . 32  0  2 
26.2435 


27.1681 

28.0940 

29.0208 

2Q.9487 

30.8778 


25.60  37 

27.5945 

29.4497 

31  .8079 

26.4354 

28.4595 

30.3447 

32.7387 

27.2697 

29.3267 

31.2410 

33.6706 

28.10  63 

30.1957 

32.1387 

34.6034 

28.9443 

31.0658 

33.0376 

79.5369 

29.7849 

31.9378 

33.9378 

36.4713 

30.62  70 

32.8115 

34.8392 

37.4066 

31.4711 

33.6865 

35.7417 

39.3424 

32.3171 

34.5627 

36.6453 

39.2791 

33.1645 

35.4405 

37.5500 

40.2165 

34.01  39 

36.3197 

3R.4$5f 

41.1545 

34.8645 

37.2003 

39.3625 

42.0933 

35.7179 

38.0817 

40.2702 

43.0326 

36.5707 

38.9644 

41.1789 

4’. 9726 

3f.4?61 

39.8*84 

42.9865 

44.9132 

TABLE  1  SHEET  X 


F 

A 

I 

PROS 

•  70  000 

•60000 

•50000 

•40000 

.30000 

L 

U 

R 

XCONF 

30.00000 

40.00000 

50.00000 

60.00000 

70.00000 

E 

S 

ALPHA 

•30000 

•40000 

0 

.35  79 

.5113 

.6935 

.9170 

1.2042 

1 

1.0974 

1.3767 

1.6787 

2.0228 

2.4399 

2 

1.9138 

2.2852 

2.6743 

3.1055 

3.6159 

3 

2.7636 

3.2114 

3.6722 

4.1755 

4.7627 

4 

3.6335 

4.1478 

4.6711 

5.2369 

5.8905 

5 

4.5171 

5.0910 

5.6703 

*.2920 

7.0056 

6 

5.4107 

6.0393 

6.6697 

7.3429 

8.1112 

7 

6.3121 

6.9914 

7.6693 

8.3898 

9.2090 

6 

7.2198 

7.9466 

8.6690 

9.4341 

10.3007 

9 

8.1328 

8.9044 

9.6688 

10.4758 

11.3874 

10 

9.0502 

9.8644 

10.6686 

11.5154 

12.4697 

11 

9.9715 

10.8262 

11.6684 

12.5533 

13.5482 

12 

10.8962 

11.7897 

12.6683 

13.5895 

14.6233 

13 

11.8236 

12.7546 

13.6682 

14.6244 

15.6955 

14 

12.75  38 

13.7208 

14.6681 

15.6581 

16.7653 

15 

13.6863 

14.6881 

15.6680 

If .6905 

17.8326 

16 

14.62  09 

15.6565 

16.6679 

17.7220 

18.8977 

17 

15.5575 

16.6258 

17.6678 

18.7526 

19.9612 

18 

16.4960 

17.5960 

18.6678 

19.7822 

21.0227 

19 

17.4359 

18.5670 

19.6677 

20.8111 

22.0825 

20 

18.3774 

19.5387 

20.6677 

21 .8394 

23.1410 

21 

19.3203 

20.5111 

21.6676 

22.8669 

24.1979 

22 

20.2645 

21.4841 

22.6676 

23.8937 

25.2537 

23 

21.21 00 

22.4577 

23.6675 

24 .9201 

26.3081 

24 

22.1566 

23.4319 

24.6675 

25.9459 

27.3615 

25 

23.1043 

24.4066 

25.6675 

2*. 9711 

28.4137 

26 

24.0529 

25.3817 

26.6674 

27.9958 

29.4651 

27 

25.0026 

26.3574 

27.6674 

2^.0202 

30.5153 

28 

25.9531 

27.3335 

28.6674 

30.0441 

31.5648 

29 

26.9045 

28. 3100 

29.6674 

31.0675 

32.6133 

30 

27.8567 

29.2869 

30.6673 

32.0905 

33.661 1 

31 

28.80  97 

30.2641 

31.6673 

•^.1132 

34.7080 

32 

29.7634 

31.2418 

32.6673 

34.1356 

35.7544 

33 

30.71 79 

32.2197 

33.6673 

3*. 1576 

36.7998 

34 

31.6729 

33.1981 

3,4.6673 

’6.1792 

37.8448 

35 

32.6287 

34.1766 

35.6672 

37.2005 

38.8889 

36 

33.5851 

35.1555 

36.6672 

x  P  •  22 1 6 

39.9326 

37 

34.54  20 

36.1347 

37.6672 

3°. 2424 

40.9756 

38 

35.4996 

37.1142 

38.6672 

4". 2629 

42.0180 

39 

36.4576 

38.0940 

39.6672 

41 .2832 

43.0600 

40 

37.41 62 

39.0739 

40.6672 

42.3031 

44.1012 

41 

38.3753 

40.0542 

41.6672 

43.3229 

45.1422 

42 

39.3349 

41.0346 

42.6671 

44.3424 

46.1826 

43 

40.2950 

42.0153 

43.6671 

45.3617 

47.2225 

44 

41.2555 

42.9963 

44.6671 

46.3807 

48.2620 

-11 

42.2165 

43.9774 

45.6671 

47.3996 

49.3009 

46 

43.1778 

44.9587 

46.6671 

40.4182 

50.3396 

47 

44.1396 

45.9402 

47.6671 

4°. 4367 

51.3778 

48 

45.1019 

46.9220 

4°. 6671 

50.4550 

52.4155 

49 

46.0644 

47.9039 

49.6671 

51 .4730 

53.4530 

50 

47.02  74 

48.8860 

50.6671 

K?.4909 

TABLE  I  SHEET  2 


-’C’  -x- 


i 


PROS 

•  20  000 

.10000 

.05000 

.02000 

*CONF 

80.00000 

90.00000 

95.00000 

90. 00000 

BETA 

0.2 

0.1 

0.05 

0.02 

1.60  99 

2.3027 

2.9959 

T .9126 

2.9945 

3.8901 

4.7448 

c . 8349 

4.27  95 

5.3229 

6.2959 

7.5171 

5.5153 

6.6812 

7.7544 

n • 08  4  6 

6.72  13 

7.9937 

9.1542 

10.5809 

7.9061 

9.2751 

10.5133 

12.0276 

9.0757 

10.5325 

11.8430 

1  T. 4370 

10.2329 

11.7713 

13.1487 

14.9170 

11.3800 

12.9948 

14.4352 

16.1736 

12.5190 

14.2061 

15.7053 

17.5100 

13.65  11 

15.4068 

16.9627 

1R.8302 

METiaaaEMi 

20.1357 

15.89  73 

17.7818 

19.4430 

21 .4284 

17.01 34 

18.9583 

20.6690 

22.7096 

18.1253 

20.1283 

21.8867 

23.9812 

19.23  34 

21.2925 

23.0972 

2C • 24  3  9 

20.3381 

22.4519 

24.3012 

26.4977 

21.4396 

23.6062 

25.4992 

97.7449 

22.5384 

24.7566 

26.6919 

9° . 9846 

23.6345 

25.9026 

27.8795 

30.2183 

24.7282 

27.0454 

29.0624 

31 .4464 

25.81 96 

28.1844 

30.2408 

39.6687 

26.9089 

29.3205 

31.4150 

37.8859 

27.9963 

30.4536 

32.5856 

36.0984 

29.0820 

31.5837 

33.7528 

36.3067 

30. 1660 

32.7113 

34.9163 

37.5109 

31.2483 

33.8367 

36.0769 

3  0.7110 

32.3290 

34.9595 

37.2345 

3°. 9075 

33.40  82 

36.0801 

38.3890 

41.1003 

34.4861 

37.1985 

39.5413 

42.2901 

35.5628 

38.3152 

40.6906 

*3.4767 

36.6382 

39.4300 

41.8379 

44.6603 

37.7123 

40.5430 

42.9826 

45.8410 

38.7854 

41.6542 

44.1253 

47.0187 

39.8575 

42.7638 

45.2659 

40.1939 

40.9284 

43.8717 

46.4042 

49.3669 

41.9983 

44.9782 

4-/. 5410 

50.5373 

43.0675 

46.0832 

48.6757 

51.7051 

44.1356 

47.1869 

#9.8086 

52.8712 

45.2027 

48.2892 

50.9400 

64.0350 

46.2692 

49.3902 

52.0696 

55.1964 

TABLE  I  SHEET  3 


[  PROS 

.01000 

.00500 

.00100 

.00050 

•00010 

.00005 

XCONF 

99.00000 

99.50000 

99.90000 

9°. 95000 

99.99000 

99.99500 

BETA 

•  o.oi 

0.005  1 

0.001 

0.0005 

0.0001 

0.00005 

4.6054 

5.2984 

6.9081 

*.6009 

9.2108 

9.9037 

6.6394 

7.4310 

9.2344 

° • 998  7 

11.7575 

12.5070 

8.4062 

9.2745 

11.2297 

12.0524 

13.9290 

14.7257 

1  0.04  60 

10.9781 

13.0631 

13.9349 

15.9143 

16.7519 

11.6048 

12.5943 

14.7943 

1  * .71 02 

1 7.7826 

18.6563 

13.1087 

14.1506 

16.4566 

17.4110 

19.5680 

20.4745 

14.5712 
16.0000 
17.40  27 
18.7835 
20.1453 


21.4901 

22.8212 

24.1397 

25.4466 

26.7434 


28.03  09 

29. 30  98 
,30.58  14 
31.8459 
33.10  32 
34.3552 
35.60  07 
36.84  18 
38.0773 

39.30  80 


40.5348 
41.7571 
42.9754 
44. 18  99 
45.40  ;'8 
4  6.60  86 

4  7.81  30 
49.0144 
50.21 28 
51.4083 
52.6010 

5  3.79  14 
54.9793 
56.16*7 
5  7.  34  78 


58.52  86 
59.7070 
60.88  38 
62.0585 
63.23  11 


64.4(1  16 
65.57  nq 
66.7382 
67.90  34 
69.0675 


15.6604 

17.1343 

18.5787 

19.9985 

21.3979 


22.  1  f'i  f 
24.1456 
25.4968 
26.8365 
28.1646 


29.4823 
30.7908 
32.0909 
33.3833 
34.6685 
35.9468 
37.2186 
38.4846 
39.7455 
41 .0004 


42.2514 
43.4969 
44.7389 
45.9762 
47.2  J95 
4  8.4  395 
49.6656 
50.8881 
52.1076 
53. 3241 
54.5376 
55. 74  0 1 
56.9558 
58.1 609 
59.3634 


60.5635 

61.7612 

62.9567 

64.1499 

65.3409 


66.5299 

67.7167 

68.9016 

70.0849 

71.2664 


18.0624 

19.6268 

21.1570 

22.6581 

24.1346 


27.0265 

28.4468 

29.8522 

31.2442 


2.6241 
33.9928 
35.3521 
36.7010 
38.0425 
39.3752 
40.7002 
42.0189 
4,3.3309 
44.6366 


45.9364 

47.2307 

48.5197 

49.8037 

51.0834 

52.3587 

53.6294 

54.8957 

5'-.  1589 

5  7.41 18_ 

58.6735 

59.9255 

61.1743 

62.4199 

63.6626 


64.9019 
6  6  •  1  3  Q 1 
67.3731 
68.6042 
69.8335 


71.05 
72.28  37 
73.5053 
74.7249 
76.9423 


19.0555 
20.6549 
22.2173 
9*  .7500 
2* .2563 


r7 .7401 
?e  .2035 
2  r>  •  6  5  0  7 
31 .081 3 
*0.4978 


**.9018 
tt .2942 
’< .6761 
*  P  •  0  4  7  9 
*r.4100 
40.7644 
42.1101 
4*  .4492 
44.7806 
4<  .1055 


4*.4248 
40.7380 
co.0456 
*1.3479 
*o.6453 
*3.9378 
**  .2258 
*4.5093 
*9.7890 
5n.  0649 
60.3368 
61 .6046 
62.8697 
64.1311 
6*.  3890 


66.6443 
67.8959 
69.1453 
70.3913 
71  .6351 


'  .8761 
74.1142 
7*. 3504 
76.5837 
77.8149 


22.9633 

24.5949 

26.1932 

27.7630 


29.3067 

30.8293 

32.3319 

33.8167 

35.2860 


36.7413 

38.1829 

39.6127 

41.0317 

42.4403 

43.8393 

45.2292 

46.6107 

47.9847 

49.3516 


50.7108 

52.0642 

53.4107 

54.7521 

56.0871 

57.4172 

58.7424 

60.0624 

61.3777 

fe2tfeflS5- 

6  3.9951 
65.2977 
66.5964 
67.891  4 
69.  1828 


70.4707 

71.7551 

73.0362 

74.3139 

76.5887 


76.8609 
78.1299 
79. ^958 
80.6598 
81.9206 


23.9298 

25.5891 

27.2135 

28.8077 


30.3762 

31.9213 

33.4464 

34.9530 

36.4431 


37.918 
39.3797 
40.8289 
42.2665 
43.6929 
45.1098 
46.5176 
47.9164 
49.3070 
50.6901 


52.0663 
53.4353 
54.7973 
56.154? 
57.5042 
68.8497 
60.1890 
61.5237 
62.8536 
64,1785 
65.4988 
66.8154 
68.1278 
69.4361 
70. 7406 


72.0414 

73.3387 

74.6325 

75.9230 

77.2102 


44 

79.7759 

81.0545 

«2.*300 

83.6024 


II 


TABLE  1  SHEET  4 


ooooo  o  c*oo 


PROGRAM  POISSN 

SC  INPUT «  OUTPUT,  TAPE5=INPUT,  T  APE6  =  0UTtrUT  ) 


COOED  F OR  RHODEN  BY  RAOKE  ON  P5MAR82. 

REAL  ANSC51,8> 

REAL  0(51) 

REAL  DELTA 
REAL  FACTOR 
REAL  LAMBDA (1000) 

REAL  LTABC1000) 

REAL  PROB ( 1 00  0*51  ) 

REAL  PTABC1000) 

REAL  XP(8> 

INTEGER  I 
INTEGER  d 
INTEGER  K 
INTEGER  L 
INTEGER  LUI 
INTEGER  LUO 
INTEGER  ML 
INTEGER  MP 
INTEGER  MR 
INTEGER  NERR 
INTEGER  NF 
INTEGER  NP 

EXTERNAL  EXP 
EXTERNAL  ITERP1 
EXTERNAL  FLOAT 

DATA  LUI/5/ 

DATA  LU0'6/ 

DATA  ML/1000/ 

DATA  MR/51/ 

DATA  MP/R/ 

C 

DELTA  =  0.1 
DO  13A  L= 1 , ML 

LAMBDACL)  =  FLOATCLMDELTA 
FACTOR  =  EXP< -LAMBDACL) ) 

D ( 1 )  =  1.0 
PROBCL.l)  =  FACTOR 
DO  132  J=2,*R 

DCJ)  =  OCJ-1) *LAM8DA(L)/FL0ATCd-l  ) 
PROB  CL ♦ J)  =  PR08(L,J-1>  ♦  FACTOR*DCJ> 
132  CONTINUE 
1 3 A  CONTINUE 


Ill  CONTINUE 


READ  (UU I *51 11 1  <XP< I > ,1=1 ,MP) 

5111  FORMAT <  8E 1 0 • 3  1 

I F  C  XPm.LT.n.O  )  GO  TO  =99 

00  121  I - 1 *  MP 

I F  <  XPCn.EQ.O.Q  )  GO  TO  122 

121  CONTINUF 
NP  =  MP 
GO  TO  123 

122  CONTINUF 
NP  =  I  -  1 

Ir<  NP.EQ.O  )  Go  Tq  999 

123  CONTINUE 

DO  142  L= 1 , ML 
J  =  ML  -  L  ♦  1 
LTA8 ( L  )  =  LAMBDA ( J ) 

142  CONTINUE 

DO  148  K=1,MR 
DO  144  L=l*Ml 
J  =  ML  -  L  ♦  1 
P  T  A  B  <  L  )  =  PRO(J,K> 

144  CONTINUF 

DO  146  L=1*NP 

CALL  ITFRP1  ( VP <L ) *  PTAB*  LTAB,  ML*  2*  ' NG ( K  »L  >  »  NFRR> 

I F (  NERR.EQ.O  )  GO  TO  145 

Nc  =  K  -  1 

UR  I TF  (LUO, 6145)  NE»R*  NF*  XP(L>*  ANSC'.L) 

6145  FORMAT(lX,12HrRROR  NUMBER, I3*5X,14HFAILDRE  *  UMBER , I  3 • 5 X, 
S  11HPR0BABILI TY,F10.5,5X,3HANS»F10.5> 

145  CONTINUE 

146  CONTINUF 
148  CONTINUF 


I 


c 

LURD  =  5 
LUWT  =  6 

y R I TF  (LUWT,fei«l) 

6181  F0RMAT<1H1*2X,1HF) 

yRITF  (  LUWT  *6  IP  2  )  <  XP  ( 1)  , 1  =  1 «  \'P  ) 

6182  FORMAT  <  3X , 1  OH  A  PRQB  •  fl  (  1  V  «  F,9  .  5  > ) 
WRITE  (LUWT, 6183) 

6183  FORMAT  <  3  X  ,1HI) 

WRITE  (LUWT.6184) 

6184  FORM AT ( 3X  ,1HL> 

DO  152  I  =  1 ,  NP 

x p < t >  =  d.o  -  xP<m*ioo.o 
15?  CONTINUE 

WRITE  (LUWT  ,6  185)  < XP <  I) «  I  =  1  * ' P  ) 

6185  FORMAT ( 3X  *  1  OH U  XCONF , 8 <1 X  ,  R9 . 5  )  ) 

WRITE  ( LUWT  *  6  186 ) 

6186  FORMAT<3X  ,1HR) 

WRITE  (LUWT, 6187) 

6187  FORMAT ( 3X  ,1H* ) 

DO  154  1=1, NP 

XP<I )  =  XP< I) /100.G 
154  CONTINUE 

WRITE  (LUUT, 6188)  ( XP ( I ) , I = 1 , XP ) 

6188  FORMAT ( 3X , 1  OH  R  ALPHA, 8( 1 *  *F9 .5  >  ) 
WRITE  (LUWT, 6193) 

6193  FORMAT (1HO) 

C 

DO  192  K= 1 , MR 
NF  =  K  -  1 
C 

WRITE  (LUO, 61^4)  NF,  ( ANS ( K , L ) , L= 1 , NP ) 

6194  F0RMAT(1X,I3,'X,8(1X,P9.4)) 

C 

192  CONTINUE 
C 

6.)  TO  111 
C 

999  CONTINUE 
STOP 
END 


subroutine;  iterpi 

*<X»XT*YT, NXtNN,Y,IER> 


C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


CODED  BY  ROBERT  RAOKE/DRSMI-RKA  ON  31MAR82 


INPUT 

X  THE  INDEPENDENT  VARIABLE  VALUE 
XT  THE  TABLE  OF  INDEPENDENT  VARIABLE  VALUES 
MUST  RE  IN  ASCENDING  ORDER 

YT  THE  TABLE  of  DEPENDENT  VARIABLE  VALUES  YT{I>  =  F  ( XT ( I ) ) 
NX  THE  NUMBER  OF  POINTS  IN  XT 

NN  THE  NUMBER  OF  POINTS  USED  IN  THE  INTEROLATING  POLYNOMIAL 
OUTPUT 

Y  THE  INTERPOLATED  DEPENDENT  VARIABLE  VALUE  Y=FtX> 

I ER  THE  ERROR  CODE  WHERE 

=  0  IF  XT  (1 )  .LE.  X  .LE.  XT(NV) 

=  1  IF  X  .LT.  XT ( 1 ) 

=  2  IF  X  .GT.  XT  C  NX ) 


REAL  XT (  NX  ) 
REAL  YT (NX ) 
REAL  X 
REAL  Y 
REAL  P 

INTEGER  L 
INTEGER  R 
INTEGER  M 
INTEGER  D 
INTEGER  NX 
INTEGER  NN 
INTEGER  I ER 
INTEGER  NP 
INTEGER  IH 
INTEGER  I 
INTEGER  J 
INTEGER  K 


noon  non  noon  non  oor»r» 


ENSURE  THAT  THE  NUMBER  OF  POINTS  IN  THE  INTERPOLATING  POLYNOM 
IS  NO  GREATER  THAN  THE  NUMBER  OF  POINTS  IN  THE  TABLE 

IF<  NX.LT.NN  )  GO  TO  112 

NP  =  NN 
GO  TO  113 

112  CONTINUE 
NP  =  NX 

113  CONTINUE 

CHECK  FOR  X  OFF  THE  LOWER  END  OF  THE  TABLE 

j F (  X.NE«XT<  1)  )  GO  TO  1 ?1 

T  =  YT<  1) 

IER  =  0 
GO  TO  999 

121  CONTINUE 

I F (  X.GT.XTC  1)  )  GO  TO  122 

X  IS  OFF  THE  LOWER  ENO 

MAKE  Y  HAVE  THE  VALUE  OF  THE  LOWER  ENO 
RAISE  THE  ERROR  FLAG 

Y  =  YT (  1) 

IER  =  1 

GO  TO  999 

CHECK  FOR  X  OFF  THE  UPPER  END  OF  THE  TABLE 

122  CONTINUE 

I F (  X.NE.XT<NX>  )  GO  TO  123 

Y  =  YT ( NX  > 

IER  =  0 

GO  TO  999 

123  CONTINUE 

I F  €  X.LT.XT(NX)  )  GO  TO  12A 

X  IS  OFF  THE  UPPER  ENO 

MAKE  Y  HAVE  THE  VALUE  OF  THE  UPPER  END 
RAISE  THE  ERROR  FLAG 

Y  =  YTCNX) 

IER  =  2 

GO  TO  999 
12A  CONTINUE 
IER  =  0 


o  n  o 


C  PERFORM  A  BINARY  SEARCH  TO  LOCATE  *  IN  XT 

C 

L  =  1 
R  =  NX 

131  CONTINUE 

M  =  (L  ♦  R>/2 

I F  <  X • NE • XT (  Y)  )  GO  TO  132 

Y  =  YT(  M) 

GO  TO  999 

132  CONTINUE 

I F (  X.LT.XTt  M)  )  GO  TO  133 

L  =  M 
GO  TO  134 

133  CONTINUE 
R  =  M 

134  CONTINUE 
D  =  R  -  L 

I F (  O.NE.l  )  60  TO  131 

C 

C  DETERMINE  THE  BEG I N I NG ( LE F T )  AND  EMn I NG ( R 1 GHT )  TABLE  POSITIONS 

C 

IN  =  NP/2 
K  =  IH  ♦  1 

I F (  K.GT.R  )  GO  TO  142 

I  =  R 
GO  TO  143 

142  CONTINUE 
I  =  K 

143  CONTINUE 

K  =  I  -  IH  ♦  NP  -  1 
I F  <  K.GT.NX  )  GO  TO  152 

L  =  I  -  IH 
R  =  L  ♦  NP  -  1 
GO  TO  154 
152  CONTINUE 

L  =  NX  -  NP  ♦  1 
R  r  NX 

154  CONTINUE 

EVALUATE  THE  INTERPOLATION  POLYNOMIAL 

Y  =  0.0 
DO  168  J=L»P 
P  =  1.0 
DO  164  1  =  Lt  R 

IF(  I.EQ.J  )  GO  TO  163 

P  =  p*(X  -  XT  <I))/(XT(J)  -  XT(I)» 

163  CONTINUE 

164  CONTINUE 

Y  =  Y  ♦  P*YT<  J> 

168  CONTINUE 

C 

999  CONTINUE 
RETURN 
ENO 


PROBABILITY  OF  FAILURE 


.07 


CONKIIJKNCK 


PROBABILITY  OF  FAILURE 


2,000 


CONFIDENCE  X 


10,000  1—  99.9  .001 


PERCENT  CONFIDENCE 


PROBABILITY  OF  FAILURE 


✓  1 


CONFIDENCE  IN  PERCENT 


HUMBER  TESTED 

PROBABILITY  OF  FAILURE  OR  TEST  TIME 


900*0 
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